Coherent reaction of Fermi superfiuid on correlated disorder 
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Motivated by the disparity between the experimentally observed properties of superfiuid 3 He in 
aerogel and predictions of the Abrikosov and Gorkov theory of superconducting alloys we consider 
effect of correlated pair-breaking impurities on the temperature dependence of the square of mag- 
nitude of the average order parameter and related thermodynamic properties of a Fermi superfiuid. 
We show that the correlations, increasing the spectral density of the fluctuations with small wave 
vectors increase the transition temperature and enhance a temperature region below T c where fluc- 
tuations of the order parameter dominate. Outside of this region we treat effect of impurities as a 
perturbation and express corrections to thermodynamic properties in terms of the structure factor 
of impurities. Assuming a simple model expression for correlation function of aerogel we find cor- 
rections to the temperature dependence of the NMR frequency shifts for A-like and B-like phases 
of 3 He and compare these with experimental data. 

PACS numbers: 67.30.hm, 67.30.er 
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I. INTRODUCTION 

Effect of impurities on traditional superconductors is 
well understood within the Abrikosov and Gor'kov (AG) 
theory of superconducting alloys^ For some of more re- 
cent superconducting and superfiuid Fermi systems, in- 
cluding high-Tc compounds and superfiuid 3 He, applica- 
tion of this theory is not always justified. Superfiuid 3 He 
in high porosity silica aerogel 3 , is a good example of such 
situation. Thermodynamic properties of this object show 
both quantitative and qualitative discrepancies with the 
generalization of the AG theory for the p-wave Cooper 
pairing^ The observed difference between the temper- 
ature of superfiuid transition in the bulk liquid and in 
aerogel T c \, — T ca at pressures above 20 bar is about a 
half of the value, predicted by the theory^ As an exam- 
ple of qualitative discrepancy we can cite the anomaly 
in the temperature dependence of the superfiuid density 
and of the magnitude of transverse NMR frequency shift. 
In Ginzburg and Landau region both quantities instead 
of the linear growth with T ca — T demonstrate the de- 
pendence which can be approximated by a power law 
(T ca — T) n with n « 1.3 — 1.4 and this dependence ex- 
tends well below the T c . 6 

The model calculations of the structure of aerogel show 
that it is an object, intermediate between the ensem- 
ble of conventional impurities and a porous mediai^ 
On the one hand the scattering centers having diame- 
ter d sa 3 nm are much smaller than the coherence length 
in 3 He (£o « 20-^80 nm), on the other - the structure has 
low density regions ( "voids" ) with the characteristic size 
of the order of £o- For improvement of the agreement 
with experiment more involved models were suggested 
which combine elements of description in terms of pores 
with that of impurities - isotropic inhomogeneous scat- 
tering model (IISM){2 and its simplified version - phe- 
nomenological IISM™, which is an interpolation between 
the two limits The models contain additional parameters 



(only one in the simplified version) . The use of these as 
fitting parameters renders a satisfactory agreement be- 
tween the calculated properties of 3 He and experimental 
data£ Still these models have certain shortcomings on 
a theoretical side. The additional parameters are intro- 
duced on an intuitive basis without precise definition of 
their relation to the structure of aerogel. As a result the 
values of these parameters can not be found indepen- 
dently. It is not clear also how accurately predictions, 
based on these models describe effect of real aerogel on 
3 He. 

In the present paper we also take into account inho- 
mogeneity of aerogel. Structural fluctuations, i.e. ran- 
dom deviations of the local density of impurities from 
the average density are treated as a perturbation. We 
concentrate on the effect of these fluctuations on ther- 
modynamic properties of the superfiuid 3 He near the T c . 
A starting approximation is the AG theory, where impu- 
rities are treated as a continuous background, producing 
damping of quasiparticles and characterized by one pa- 
rameter - mean free path of quasiparticles I. The first 
correction to the AG approximation can be expressed in 
terms of the structure factor of aerogel, which has precise 
definition and can be measured directly. The obtained re- 
sults arc approximate, but their accuracy is controlled by 
the small parameter, which is a combination of parame- 
ters, characterizing both aerogel and superfiuid 3 He. 

For account of structural fluctuations of pair break- 
ing impurities we follow the argument of Larkin and 
Ovchinnikovi^i They have shown that on a low temper- 
ature side of the transition (ip) 2 in addition to the term 
depending linearly on (T c — T)/T c acquires a singular cor- 
rection proportional to y/(T c - T)/T c . At T -> T c this 
correction can catch up with the principal term. That 
renders an estimation of a region of "broadening" of the 
transition ST{,r- At (T c — T) <C ST(, r temperature de- 
pendence of p s and other thermodynamic properties of 
a superfiuid is determined by fluctuations For Gaus- 
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sian fluctuations 6Tb r ~ (T c /N 2 )(^q/1) 4 , where N% is a 
number of scattering centers within a volume with the 
radius ~ £q- In aerogel the cross-section of a scatter- 
ing center is ~ 10 2 times greater than that of an atomic 
impurity in a metal. It means that for the same mean 
free path number of centers in a unit volume is ~ 10 2 
smaller. Coherence length £o for 3 He is of the same or- 
der as in a typical superconductor, so the fluctuations 
are bigger. Still for a 98% aerogel N% ~ 50 and con- 
tribution of structural fluctuations dominates only in a 
region (T c - T)/T c ~ 10 ~ 4 - 1CT 5 , which is too narrow 
to explain the observed anomaly. 

Porto and Parpia 7 pointed to the fact that in aerogel 
the distribution of scattering centers is far from being 
random. They gave convincing experimental evidences of 
strong effect of structural correlations within the aerogel 
on the shift of the T c and on the temperature dependence 
of p s . Guided by these observations we consider theoreti- 
cally effect of correlations in distribution of pair-breaking 
impurities on the shift of the T c and on the temperature 
dependence of the square of the absolute value of the or- 
der parameter of a Fermi superfluid in a vicinity of T c . 
We introduce in the argument of Ref. [ll| effect of corre- 
lations. Because of the coherence of condensate its inter- 
action with the ensemble of impurities depends on phase 
shifts between the perturbations brought about by the 
centers. The result of interference brings in the numera- 
tor of the expression for broadening of the transition an 
extra factor ~ (nR 3 ) 2 , where R is the correlation radius. 
At sufficiently large R this factor compensates for N? in 
the denominator and effect of structural fluctuations can 
extend over an appreciable region below the T c altering 
in this region temperature dependence of thermodynamic 
properties of the superfluid. 

The paper is organized as follows: in Section [TT] we 
consider correlated impurities in a superconductor with 
a scalar order parameter. Short account of these results 
was published before i 13 ' 14 Here we added as illustration 
an example of a system with a realistic correlation func- 
tion and dropped all results of Ref. [U, which go be- 
yond the perturbation theory and contain uncontrolled 
approximations. Reformulation of the argument for the 
p-wave superfluid 3 He is presented in Section HTT1 In Sec- 
tion IIVI comparison with the NMR data for ABM and 
BW phases of 3 He is made. 



II. SCALAR ORDER PARAMETER 

In a vicinity of T c effect of a quenched disorder can 
formally be described as a random spatial variation of 
coefficients in the Ginzburg and Landau equation. For a 
scalar order parameter \E': 



o(r)$ + 6(r)*|*| 2 - V(c(r)W) = 0. 



(1) 



Following the procedure of Ref. |Tl| we take into account 
only spatial dependence of a(r), because it renders the 



most singular contribution to the thermodynamic prop- 
erties at T — > T c . For 6(r) and c(r) we use their av- 
erage values (b) and (c). The coefficient a(r) can be 
rewritten in terms of the local transition temperature 
T c (r): a(r) = a(T — T c (r)). It is convenient to divide 
T c (r) into its average value (T c ) and a relative fluctua- 
tion r/(r): a(r) = a(T c )[r- T](r)], where r = -1+T/(T C ) 
and t)(t) = -1 + T c (r)/(T C ). After substitution * = * V-> 
with = a(T c )/(b), and £ 2 = (c)/a(T c ) Eq. (1) takes 
the form 



[r - t]( T )]i/> + 4>\4>\ 2 - £ 2 A^ = 0. 



(2) 



The random function 77 (r) is treated as a perturbation. 
The global T c is defined by the condition (ip) 7^ at 
T < T c . Below T c solution of Eq. (2) can be sought 
in a form ip( r ) = (V'Hl + x( r )) an d the perturbation 
procedure is justified when (x 2 (r)) <C 1. For real rj(r) 
essential part of \ is also real. Keeping in Eq. (2) terms 
up to the second order in x and r\ we arrive at: 

t- v + t X -VX+ W 2 (l + 3x + 3 X 2 ) - C 2 Ax = 0. (3) 

By the definition (r](r))=0 and (x(r))=0. Taking average 
of Eq. (3) we arrive at the equation for (ip) 2 : 



w 2 



(vx) - 1 

l + 3(x 2 )' 



(4) 



Subtraction of the averaged part from Eq. (3) renders 
linear equation for x( r ) : 

(r - 77(F) + 3W 2 )x - £ S 2 A X = r?(r) - ( VX ). (5) 

The second order term 3(-0) 2 (x 2 — (x 2 )) is neglected in 
comparison with 3(ip} 2 x- Equation © can be formally 
rewritten in terms of the Green function G(r, r'): 



X(r)= / G(r,r')(r/(r')-(ryx))rf 3 r'- 



(6) 



A straightforward argument shows that the average (rjx) 
can be expressed via the self-energy of the averaged 
Green function of Eq. ([5]) in the momentum represen- 
tation S(k, t): 



(XT?)=E(0,r). 
In a principal order on the perturbation ?y(r) 



£(0,t) 



M-kMk)) 



d 3 k 



+ 3(^) 2 -S(k,r)+C 2 fc 2 (2ir) 3 ' 



(J) 



(8) 



Now we take u — S(0, r) — r as a new variable. The 
self energy in the denominator of Green function may be 
expanded as E(k, t) sa S(0, r)+ 3 ^ fc 2 and may be 

included in the definition of £ 2 Transition temperature is 
determined by the condition u = 0. In the second order 
on rj(r) it is 



Tr>. 



(T c )[l 



fa(-k)Tj(k)) d 3 k 



gk 2 



(2tt) 3 



(9) 
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Using these definitions and Eq. (8) we arrive at the rela- 
tion between u and t=(T— T c2 )/(T C ): 



1 



(r?(-k)r?(k)) 



d 3 k 



e s V[2u + Qk 2 ] (2tt) 



-t. 



(10) 



In evaluation of the average (?7(-k)?7(k)) we neglect effect 
of multiple scattering of quasiparticles by the scattering 
centers. That introduces a relative error in evaluation of 
(T c ) of the order of x 2 / 3 , where x is the volume concentra- 
tion of impurities. Assuming also that all scattering cen- 
ters are identical we have: r)(r) = *7^( r— r a)~ if 1 (0) 
, where rf 1 '^ — r a ) is a local suppression of transition 
temperature by one center situated at the position r Q , 
7y( 1 -'(k) - its Fourier transform , t/^O) = //^(k = 0). 
Then (77(-k)7y(k)) = n\r]^ (k)| 2 S'(k), where the structure 
factor S'(k) = exp[ik(r;, — r a )]) depends on the dis- 
tribution of impurities. If impurities are not correlated 
all terms with ^ r a for finite k vanish at the averaging. 
Then S'(k) = 1 + (27r) 3 n£(k). Here unity comes from the 
summand with r& = r a . The term proportional to 5(k) 
can be dropped because it eventually enters with the fac- 
tor (r](r)) = 0. When correlations are present terms with 
7^ r a render finite contribution to the structure factor. 
The number of additional terms in the sum is of the or- 
der of nR 3 , where R is correlation radius. If nR 3 ^> 1 
their sum, which in what follows is denoted as S*(k) can 
be a principal contribution to S(k). The S'(k) is directly 
related to the correlation function in the coordinate rep- 
resentation. For globally isotropic distributions the prob- 
ability w(rb\r a ) to find a particle in the point rf, if there 
is a particle in the point v a depends only on a distance 
r = |r;, — r |. At r — > oo correlations vanish and w(r) 
tends to a constant. Normalization of w(r) is usually 
chosen so that this constant is unity. Then a measure of 
correlations is v(r) — w(r) — l. Changing in the definition 
of S(k) summation for integration we arrive at 



S(k) 



v(r)e' 



- lkr d 3 r 



Eq. (fit))) can now be rewritten as 

u[l + W{u)\ = -t, 



where 



W{u) 



n|? ? W(k)| 2 [l + S(k)] d 3 k 
e s k 2 [2u + gk 2 ] (2tt) 3 



(11) 



(12) 



(13) 



Numerator of the expression under the integral is the 
spectral density of fluctuations. It is multiplied by the 
response function, which is singular at u — > and k — > 0. 
Contribution of fluctuations with k ~ l/£(u), where 
£(w) = £ S /V2u ^> £ s is strongly enhanced. The perturba- 
tion 77 W (r) falls out on a distance <~ £ s from the center so 
that in the essential region of k (k) ss n^'(0) ~ £ s er, 
where a is a cross-section of scattering of a quasi-particle 
by one center. When fluctuations are not correlated 



S(k) = 0, the spectral density is constant and we recover 
the result of Larkin and Ovchinnikov 11 for broadening of 
the transition: Wi{u) = n\v (1 \0)\ 2 /(2tt^ 3 V2u)- Rough 

For 3 He in the 
Contribution of 



estimation renders Wi(u) ~ ( y- Y 

98% aerogel W x {u) < 1 if u > 1CT 4 . 
correlations 



W 2 (u) = 2tt,|77 (1) (0)| s 



S(k) 



d 3 k 



e s k 2 [2u + e s k 2 ] (2tt) 3 



(14) 



depends on the behavior of S(k) at k ~ l/£(u). X-rays 
scattering dat a 7 ' 8 ' 15 show that silica aerogels have frac- 
tal structure with the fractal dimension Df rj 1.7 -j- 1.9 
depending on a sample. Formally it means that at k — > 
the 5(k) grows as l/kf increasing density of the long 
wavelength fluctuations until it saturates at k ~ l/R. 
Further enhancement comes from the response function. 
It is characterized by the temperature dependent coher- 
ence length £(«)• If R 3> £s there are two regions of u, 
corresponding to different asymptotic behavior of W 2 (u) : 
£(u) > R and < R with different asymptotic. For 
finding the asymptotic we express W 2 (u) in terms of the 
correlation function v(r): 

W 2 (u) = lV c2 Wl / [l-exp(-r/£(«))]i;(r)rdr. (15) 
CI" J 

In a simplest case v(r) can be characterized by two 
parameters - a correlation radius R and an overall am- 
plitude A. Convenient model form was suggested in Ref. 



v(r) 



r 



3-D, 



exp(— r/R). 



(16) 



This expression reproduces qualitative features of corre- 
lations in silica aerogels i.e. the existence of a macro- 
scopic correlation radius R and a fractal behavior with 
dimensionality Df at r < R. Introduction of the dimen- 
sionless variable y — r/R renders: 



W 2 {u) 



n 2 i? 2 |r ? ( 1 )(0)| : 



exp(-yR/£(u))]v(y)ydy. 

(17) 

Asymptotic of this integral at (R/((u)) — > is W 2 (u) ~ 
(n 2 \r]^(Q)\ 2 R 3 V2/^ 3 ) v(y)y 2 dy. It has the same 
1/s/u singularity at u — > as the contribution of non- 
correlated impurities. Its sign and the relative weight 
depend on the value of the integral = J °° v(y)y 2 dy ~ 
A. It means that the contribution of correlations has an 
extra coefficient AnR 3 . For realistic values of parameters 
effect of fluctuations in aerogel is enhanced. 

In the opposite limit R 3> £(u) the exponent in the 
square bracket under the integral in Eq. (|15|) can be ne- 
glected and W(u) v(y)ydy. That renders a 
constant shift in the relation between u and t which is 
equal to (T c2 — (T C ))/(T C ) so that eventually u « — r. 
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T/(T C ) 

FIG. 1. Dependence u on T/(T C ), as given by Eqns. (12) - 
(16) for different ratios of R/£,s at a fixed value of the product 



A{R/p) 



3-D, 



Solid line corresponds to non-correlated case, 



i.e. R/£s = 0. Dotted line corresponds to R/£ 3 
dashed line is for R/£ s = 7. 



3 and 



at 



This is a straight line, which extrapolates to T = (T ( 
u = 0. 

Relation (|T2l) between u and t was derived with the use 
of perturbation theory, it breaks down when W(u) ~ 1. 
This condition determines a region of broadening of the 
transition by structural disorder. Using the asymptotic of 
W(u) at u — y we arrive at the estimation of the width of 
this region: u < (AR 3 / £ s l 2 ) 2 . The estimation coincides 
with one obtained from the condition (x 2 ) <C 1. It follows 
from the asymptotic at R <C £ (it) of the expression 



(x 2 > = 



n|r?W(k)| 2 [l + 5(k)] d 3 /fc 



[2« + £ 2 fc 



212 



(2tt) 3 



(18) 



The region of broadening of the transition for correlated 
impurities is (AnR?) 2 times more wide than for non- 
correlated and the "tail" of the singular correction to the 
linear dependence of u on t extends for a larger interval 
of u into the region (AR 3 /£ S Z 2 ) 2 <C u <C 1 where per- 
turbation theory holds. If R ~ £ s the deviations of (ip) 2 
from the linear law can persist in a region well below the 
transition temperature in agreement with the data for 
superfluid 3 He in aerogel^ 

For illustration of qualitative changes, introduced by 
correlations we plot on the Fig. [T] the dependence 
u{T/{T c )), given by Eqns. (12) - (16) for different val- 
ues of parameters, characterizing aerogel. As one can 
see the transition temperature shifts up with respect to 
(T c ) when R increases. Relative value of the shift ac- 
cording to Eq. ([9]) is of the order of (R/l) 2 . Ratio of the 
positive shift to T c b — (T c ) is of the order of R 2 /^ s l. If 
R ^> £ s a. major change of the slope of u(t) takes place 
at u ~ (6/i?) 2 . 

Function u(t) determines also correction to the temper- 
ature dependence of the specific heat below T ca . Differ- 
ence of free energies of the superfluid and normal phases 



can be represented as: 



F s -F n = -^(T C ) 2 



2T, 



cb 



(19) 



where ACj is the specific heat jump in pure 3 He. With 
account of the terms of the order of (\) 2 (|V'| 4 ) = u2 > 
then: 



T 



(C s - C n ) 



AC b d 2 (u 2 ) 
2T cb dt 2 



(20) 



For u{t) given by Eqns. (12), (13) this expression remains 
finite up to T ca , but its validity is justified only in a region 
where the condition (% 2 ) 1 is met. 



III. P-WAVE PAIRING 

For application to the superfluid 3 He in aerogel the 
argument of Section |TT] has to be reformulated for the p- 
wave Cooper pairing. The order parameter in that case is 
3x3 complex matrix A^j . The first index refers to 3 spin 
components, and the second - to 3 components of the an- 
gular momentum of Cooper pairs. Aerogel interacts with 
the orbital part of the order parameter. This interaction 
is described by the additional term in the density of free 
energy: 



fv 



(21) 



Now rjji(v) is a real and symmetric 3x3 random ten- 
sor of local anisotropy. The ensemble average (rjji) = 0. 
The r-independent tensor kji is introduced for account 
of possible global anisotropy. In the present calculations 
a uniaxial anisotropy will be assumed only implicitly for 
lifting of the orientational degeneracy of the A-like phase. 
Within the used approximation the global anisotropy will 
not enter results and it will be dropped from the free 
energy. N e ff is the overall factor, having dimensional- 
ity of density of states, it does not enter the resulting 
equations. For the gradient energy we take a simplified 
isotropic expression: 



/v = N eff i 



2 ( djna 9A^ 

V dx n 8x n 



(22) 



Both terms (21) and (22) are added to the unperturbed 
density of GL free energy: 



fo = NeffirA^A^ 



(23) 



s=l 



where -L 



are the 4-th order invariants in the expan- 

.16 



sion of the free energy over A^jr^ I\ 



— A,,i A 



M^ljAviA*^, 



A , A* 



— A A ■ A* A* 



A* 



L 5 — ■"■lij'-uj 

phenomenological coefficients, 
tions have the form: 



h = A.jAlA^A*,, and (3 



, h = 
.../3 5 are 



The resulting GL equa- 



dfo 
dA 



1 1 3 



d 2 A 



dxl 



= -A^ mj (r). 



(24) 
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To follow the perturbation procedure of the previous 
section wc split the order parameter into its average and 
the fluctuation. In the B-like phase we search for a 
solution in the form A^j = {A^i)(5ij + cy(r)). Here 
{Aftj} = Ai? M: j, R^j is a real orthogonal matrix and cy(r) 
is real and symmetric. Procedure analogous to that used 
at the derivation of Eq. (0| renders an equation for A 2 : 

2 1 

+ 12 (c mm cii)} = -(rimicim) - T. (25) 

Here /3i2 = Pi+02 etc.. Tensor of fluctuations Cji satisfies 
the linear equation: 

(rSjn - r) jn )c nl - £ 2 V 2 c 3 i + A 2 [3/3 12345c 3 i 

+/3i 2 2c nn Sji] = Tfji - (rjjnCni). (26) 

As before we consider aerogel as a random distribution 
of identical scattering centers. The center situated at a 
point r Q introduces a local perturbation m]' (r — r ). For 
sufficiently dilute system, when effect of multiple scatter- 
ing of quasi-particles by impurities can be neglected, the 
overall perturbation is a sum of contributions of indi- 
vidual impurities: r]ji(r) = ^2 a Vji ( r ~ r a)- Its Fourier 
transform is: 

^(k) = r ? «(k)^expHkr a ). (27) 

a 

As has been discussed in Section [H] most essential are 
fluctuations with small k, which meet the condition 

fc£ s "C 1. That makes possible to substitute for rf^ (k) 
its limiting value at k — > which is c^r/ 1 ' (0), leaving 
dependence on k only in the phase factors exp (— ikr a ). 
In the coordinate representation it means that the scat- 
tering centers are considered as isotropic and point-like 
r]jl\r) — 5jir)^(0)S(r — r a ) The total tensor 77^7 (r) in 
that case is also proportional to Sji. By the argument 
of Rainer and Vuorioii r]^(0) = 7(£ a /n£) where I is a 
transport mean free path and the coefficient 7 depends 
on properties of the centers. For diffusely scattering balls 
7 = 7r 2 /4ji2 Symmetric tensor Cji in Eq. (f26| can be 
represented as a sum of its scalar c nn 5ji and traceless 

c ^ji = \i c J l + Cl i ~ ^ c nn5ji) parts. The system of equa- 
tions (|2"6"|) splits then in two independent systems for each 
part of Cji. Of the two only equation for the scalar part 
contains a finite perturbation in the right hand side, so 
that the solution can be sought in a form Cji = SjiX- Sub- 
stitution of rjji = r]6ji in Eqns. (25), (|26p and introduc- 
tion of the notation (ip) 2 = A 2 (3/?i2 + transforms 
Eq. (25) into Eq. Q and Eq. flSBJ) - into Eq. ©. The 
problem for 3 He-B turns out to be equivalent to that for 
the scalar order parameter. 

In the A-like phase situation is complicated by a pos- 
sibility of formation of the disordered Larkin-Imry-Ma 
(LIM) state. Global anisotropy orients the order param- 
eter and restores the long-range order. Experiments 1 ^ 
and theoretical estimation a 19 ' 20 show that compression 



of aerogel for a few percent along one direction is suffi- 
cient for the restoration of the long range order. In this 
case reduction to the scalar problem is possible for the 
A-like phase as well. In the bulk liquid the order param- 
eter of the A-phase is a direct product of spin and orbital 
vectors: A^j — d^Aj. Aerogel coated by 4 He does not 
interact with the spin part of the order parameter. We 
seek a solution of Eq. (24) also in a separable form, then 
it reduces to the equation for the orbital vector Aj: 

[rSjt - Vj i(r)]Ai e s (^A + PisA*(A s A s ) 

+0 MS A j {A.Al)=O (28) 

In line with the discussion in the paragraph follow- 
ing Eq. (27) we keep only isotropic part of the pertur- 
bation r)ji(r) = Sjir}(r) and seek the solution in a form 
Ai = (Ai) + ai. Here (A t ) = (l/v / 2)A(m + in), unit 
vectors rh and h are mutually orthogonal and orthogonal 
to the direction of compression. For definiteness we have 
chosen one of two possible orientations of (A{). Repeat- 
ing the procedure of Section |TT] we find as the solution 
a; = x(r)(A;). Denoting /3245A 2 = (ip) 2 we return to the 
equations (4) and (5) of the scalar problem. 

The reduction to the scalar order parameter is based 
on the fact that principal contribution to thermody- 
namic anomalies comes from inhomogeneities with small 
wave vectors. Such inhomogeneities interact with fluc- 
tuations of the amplitude of the order parameter, their 
coupling with the other collective modes of B-like and 
A-like phases appears only in the next order on u. 



IV. COMPARISON WITH EXPERIMENT 

The obtained expressions allow comparison with ex- 
periment. Properties of aerogel enter these expressions 
via the mean free path / and the structure factor S'(k). 
Both quantities in principle can be found from indepen- 
dent experiments, I - from the measurement of spin dif- 
fusion in the normal phased and S(k) - from the X- 
ray scattering data. Physical quantities which can be 
compared with our calculations are e.g. the temperature 
dependence of the square of the longitudinal resonance 
frequency f2 2 , which in a vicinity of T c is proportional to 
A 2 i.e. (^) 2 (Eqns. (4), (12), (13)) and the lowering of the 
transition temperature, given by Eq. ([9]). In realization 
of the comparison we encounter two difficulties. Firstly 
the thermodynamic data suitable for the comparison are 
not accompanied by the independently measured I and 
S'(k) for the same sample of aerogel. Even when 5(k) is 
measured the result is given in arbitrary units, leaving 
ambiguity in the overall amplitude. Secondly the theo- 
retical expressions are obtained within the perturbation 
theory and their application is limited by a region of small 
fluctuations. Since the discussed effects originate from 
fluctuations they are small in a region of applicability of 
theory. So, we have to find a compromise between the 
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FIG. 2. Fitting of the dependence of lowering of the transition 
temperature of 3 He in aerogel (data points from Ref. 0) as a 
function of £o = 2-kHvf /T cb by the straight line "f£,o/l + 8 (solid 
line). The slope 7 = 2.95 falls between its values for specular 
7s P = 7r 2 /3 and diffuse scattering 7^/ = 7r 2 /4. The constant 
shift 8 — 0.075 is a contribution of correlations to the change 
of the transition temperature (8 ~ AR 2 /I 2 ). The dotted line 
corresponds to non-correlated impurities. 



magnitude of the effect and accuracy of its description. 
With these reservations we choose for comparison a tem- 
perature dependence of fl 2 ^, extracted from the measure- 
ments of the transverse NMR frequency shift. 22 In these 
experiments neither I nor S(k) was measured, but the 
data was taken for several pressures with the same sam- 
ple of 98% aerogel and in a suitable temperature interval. 
For such porosity I is roughly estimated as 130-180 nm. 21 
Fine tuning of I within this interval is made in a process 
of fitting. Unknown structure factor formally introduces 
infinitely many fitting parameters, but in all equations 
5(k) enters under the sign of integral. The main contri- 
bution to the integral in Eq. (13) comes from the region 
of small k so that only the asymptotic form of 5*(k) at 
k — > or of v(r) at r — > 00 is essential. The integral in 
Eq. © for Df > 1 is also defined by a region of small k. 
The answers are not very sensitive to the detailed form 
of S(k) or of v(r), so we use for interpretation of data 
the model expression Eq. (16) with two fitting parame- 
ters A and R. Lowering of the transition temperature 
by aerogel was discussed before*^ In the leading order 
on £oA the principal contribution is proportional to this 
ratio and correction due to correlations depends on the 
combination AR 2 /l 2 . For comparison we use the data of 
Ref. 11 taken at four pressures: 29.3, 17.5, 20.1 and 11.9 
bar. Transition temperatures for all four pressures can 
be fitted by a straight line as shown on the Fig. [2] and 
the combination AR 2 /I 2 can be extracted from the fit. 
Additional combinations of parameters are found by fit- 
ting temperature dependencies of f2 2 in the B-like phase, 
as shown on the Fig. [3] The procedure requires further 
fitting parameters originating from the phenomenologi- 
cal coefficients, entering Eqns. (1), (24). The amplitude 
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T/T ca 

FIG. 3. Temperature dependence of the square of the longi- 
tudinal resonance frequency Q% in the B-like phase at three 
different pressures. Data points are taken from Ref. [f| solid 
lines - results of fitting. Parameters of aerogel are I — 140 nm, 
R = 35 nm, A = 0.055. Coherence length is taken as 
£ s = 0.71&. 



has to be fitted for each pressure. About the coher- 
ence length £ s we assume that it scales as {hvF)/T c b and 
fit one coefficient for all pressures. Parameters of aero- 
gel, which provide the best fit for the lowering of the 
transition temperature for all four pressures and for the 
temperature dependence of the square of the longitudi- 
nal resonance frequency £l 2 L for three pressures, shown 
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FIG. 4. Temperature dependence of the square of the longitu- 
dinal resonance frequency Q,\ in the A-likc phase at P=29.3 
bar. Data points are taken from Ref. |23 . solid line - results 
of fitting with parameters of aerogel - I = 140 nm, R = 35 
nm, A — 0.055. Coherence length is taken as £ 3 = 0.71£o. 

on the Fig. [5] are the following: I — 140 nm, R = 35 nm, 
A = 0.055. 

Part of the data fall in a region were (x) 2 is not 
very small. Estimation, made with the cited parame- 
ters renders for P=11.4 bar (x) 2 » 0.04/ y/u at u -> 0. 
At the position of the closest to T c experimental point 
(x) 2 ~ 1/3- Fluctuations are even bigger at P=20.1 
bar, where (x) 2 ~ 0A/y/u. Argument, based on the 
perturbation theory can not apply in a region of devel- 
oped fluctuations. The argument certainly breaks down 
when (x) 2 ~ 1. For P=20.1 bar it would correspond to 
u = 0.01. 

For the A-like phase we also use the NMR data of 
Ref. [221 . taken in a condition when compression is suffi- 
ciently strong for orientation of the average order param- 
eter. Fig. |4] demonstrates that a satisfactory agreement 
with experiment can be achieved. 

The results of comparison show that the data for both 
A-like and B-like phases can be fitted with the values of 
parameters which are within the range admitted by other 
experiments with different samples of 98% aerogel. That 
provides a support to the suggested interpretation of the 



observed anomalies of thermodynamic properties of the 
superfluid 3 He in aerogel. 

V. CONCLUSION 

When discussing reaction of a superconductor or a su- 
perfluid on a quenched disorder it is important to take 
into consideration that reaction of the condensate of 
Cooper pairs increases with a decrease of a wave vector of 
a perturbation k, saturating at k ~ l/£(w). Correlations 
between impurities influence distribution of fluctuations 
over their Fourier spectra. If, like in aerogel, correlations 
increase the weight of the long wavelength fluctuations 
effect of the disorder can be significantly increased in 
comparison with the uniform spectral distribution of the 
same impurities. This results in the increase of a tem- 
perature interval below T c where fluctuations dominate. 
The fluctuation region has a "tail" , in which fluctuations 
of the order parameter are sufficiently small to be treated 
by the perturbation theory, but their effect on thermody- 
namic properties of the superfluid is still distinguishable. 
A length of the tail depends on the relation between the 
coherence length of the condensate £0 and correlation ra- 
dius R. For R £0 a trace of fluctuations extends for 
the most of the Ginzburg and Landau region. 

Shift of fluctuations in a region of small k enhances the 
transition temperature of 3 He (cf. Eq. ©). Analogous 
effect was discussed before as a way to further enhance 
the T c in high temperature superconductors by order- 
ing of impurities or modulation of paring interaction^^ 
Advantage of 3 He in aerogel as an object of investigation 
is that aerogel forms a rigid structure, which can be stud- 
ied separately from 3 He. 
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